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Abstract: The non-linear back reaction of accelerated cosmic rays at the shock fronts, leads to the for- 
mation of a smooth precursor with a length scale corresponding to the diffusive scale of the energetic 
particles. Past works claimed that shocklets could be created in the precursor region of a specific shock 
width, which might energize few thermal particles to sufficient acceleration and furthermore this precur- 
sor region may act as confining large angle scatterer for very high energy cosmic rays. On the other hand, 
it has been shown that the smoothing of the shock front could lower the acceleration efficiency. These 
controversies motivated us to investigate numerically by Monte Carlo simulations the particle accelera- 
tion efficiency in oblique modified shocks. The results show flatter spectra compared to the spectra of the 
pressumed sharp discontinuity shock fronts. The findings are in accordance with theoretical predictions, 
since the scattering inside the precursor confines high energy particles to further scattering, resulting in 
higher energies making the whole acceleration process more efficient. 



Introduction 

It has been shown theoretically 30 years ago that 
the first order Fermi acceleration, namely also dif- 
fusive shock acceleration (e.g. [1], [2]), is a plausi- 
ble mechanism for cosmic rays to gain statistically 
an amount of energy by crossing a shock front - 
formed in a super-Alfvenic plasma flow- in con- 
secutive cycles, while scattering off the irregular- 
ities of the magnetic field present in the media. 
The energy distribution of these particles follows 
a power law, which generally is in accordance with 
the actual measurable value of the spectral index 
of the observed cosmic ray spectrum, at the top of 
the atmosphere. Over the years, in the studies of 
diffusive shock acceleration it is commonly con- 
sidered that a shock front is an immaterial non- 
dimensional surface, although this is not true. It 
quickly became evident to the researchers that the 
test-particle approximation, where the particles do 
not interact with the shock, could not entirely hold 
and non-linear effects could take the leading part 
into the dynamics of the shock acceleration mech- 
anism. In other words, if a reasonable amount of 



energy was transfered to the accelerated particles 
(i.e. > 10%), they could dynamically play an im- 
portant role in the shock process itself by interact- 
ing with the shock eventually modifying it. More- 
over, the modified shock which will consequently 
have a finite width (a precursor), will have a length 
scale which will correspond to the diffusive length 
scale of the energetic particles. 
Specifically, work of [3], [12], [13] and [14] ques- 
tioned the simplification of the test-particle ap- 
proximation, which consequently considered the 
shock as an immaterial surface into the accelera- 
tion mechanism process and moreover showed that 
cosmic ray influenced shock structures are unsta- 
ble, indicating that the turbulence present could 
be strong, requiring a full treatment of a violently 
scattering medium as a result. This instability 
present in the modified shock structures, could be 
of a great importance and a numerical investigation 
following the above theoretical works, could throw 
further light into the exact mechanism, particularly 
for oblique shocks. Concerning the latter, work of 
[10] has shown that highly non-relativistic oblique 
shocks could be very efficient cosmic ray accelera- 
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tors under given circumstances of the particles dif- 
fusion [6J. Therefore, the results presented here 
could have further implications on the efficiency 
of these shocks, provided there is no immaterial 
shock surface as assumed before, but rather the 
conditions which are described in this work. Sum- 
marizing, in this work we aim to test and com- 
pare past theoretical claims with numerical simu- 
lations and investigate changes in the acceleration 
efficiency, spectral indexes and acceleration rates, 
simulating an oblique modified shock (with a fi- 
nite width). In a simplified manner, altering the 
diffusion conditions inside the precursor, this wiU 
physically reflect on the magnetic field irregulari- 
ties within the modified shock region. 

Discussion 

It is known that concerning the diffusive shock ac- 
celeration mechanism, if one wishes to treat more 
reahstically the shock acceleration mechanism it- 
self, one should take into account that the shock 
is not an immaterial surface anymore and the text 
particle approximation can not hold. The issue 
of the back-reaction of the accelerated particles 
on the shock and the violation of the test-particle 
approximation, occurs when the acceleration pro- 
cess becomes sufficiently efficient to generate pres- 
sure of the accelerated particles, which is com- 
parable with the incoming gas kinetic pressure. 
The spectrum of the cosmic rays and the struc- 
ture of the shock are changed by this effect which 
is intrinsically nonlinear. As it is mentioned in 
the classic review of [7], one could classify the 
non-linear phenomena in different processes which 
could act on the shock region, considering that the 
acceleration is efficient. One process, is the one 
where the accelerated particles can interact with 
the unshocked plasma creating instabilities, such 
as different types of plasma waves, [11]. Another 
process, is the one where the particles built up 
a sufficient pressure P = 47r/3 / p^Vf{p,t)dp, 
which could slow the background plasma before 
the shock and modify it to an extended structure. 
One can think that as the diffusion coefficient of 
the particle is increasing with energy, then the very 
high energy particles will escape (having large dif- 
fusion lengths) the acceleration region, with a con- 
siderable energy and momentum. This in turn will 



allow for the compression ratio to increase, which 
will be followed by an increase in the accelera- 
tion efficiency. This condition will flatten the flnal 
spectrum and it wiU give the opportunity to more 
particles to get accelerated to higher energies, be- 
fore leaving the system. Furthermore, if relativis- 
tic particles carry with them most of the energy 
and produce a signiflcant part of the pressure in 
the shock system, then the shock will become cos- 
mic ray dominated and the index of the speciflc 
heats wiU decrease due to the relativistic particle 
speeds. This fact wiU initiate a further increase of 
the compression ratio, making the shock even more 
efflcient, pushing the energy towards the very high 
energy relativistic particles. Nevertheless, an in- 
crease in the compression ratio is not as easy, as [5] 
has shown with simulations and for non-relativistic 
shocks the compression ratio tends to a value of 
four, being a reasonable choice for studies even in 
the full non-linear case. 

Numerical method - Results 

In our present simplified approach, this non-linear 
problem will be dealt with in a linear way. The 
Monte Carlo numerical approach simulates as real- 
istically as possible the random process of a phys- 
ical process; namely in our case, the random walk, 
the diffusion and the scattering of particles, with 
the supposedly magnetic field irregularities, which 
can be prescribed within the simulation reflecting 
directly by the different diffusion coefficients. A 
specific feature of the simulations is (following [3]) 
that there is scattering inside the precursor (differ- 
ent from what is upstream and downstream of the 
shock formation) of length of L = k/V, where k 
is the spatial diffusion coefficient given by: k = 
K\\cos'^tp + K^sin^ij) and V is the velocity of the 
shock. The angle ip is the inclination between the 
magnetic field lines and the normal of the precur- 
sor For parallel shocks one will have L = k\\/V, 
while for perpendicular ones, L = k^/V, while 
for aU other cases it is L = k/V. In this work 
we simulate shocks of a high inclination of 75 de- 
grees (as shown in figures 1 and 2), other high in- 
clinations are comparable, since as we presented 
in [10], even for the test-particle approximation, 
highly oblique shocks can be very efficient cosmic 
ray accelerators, achieving energies of ~ 10^'^ eV. 
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Figure 1: The shock velocity versus the obtained 
spectral index for three different simulation runs 
for a shock inclination equal to 75 degrees. The 
measurements are taken in the downstream side of 
the shock rest frame. Starting from the top, the 
first Une is for a shock with no precursor, the sec- 
ond is for a shock with a precursor width equal to 
10 • A and the third one of precursor width equal 
to 100 • A. One can see that the modified shocks 
create a flatter spectra compared to the spectrum of 
sharp shock discontinuity. Also, one may see that 
the spectral index flattens as the shock's velocity 
increases. We get similar results for other inclina- 
tion angles and as long as the shocks are obUque. 

For a complete picture, in a forthcoming detailed 
paper we present a variety of shock inclinations, 
ranging from nearly parallel to perpendicular ones. 
Furthermore, since the cosmic ray diffusion tensor 
is diagonal in view of the adapted magnetic field 
orientation, one may say that in a Cartesian system 
xyz, Kxx = '^"11 'ir^d '^■yy = ^^zz = i^L, with K||, kx 
the cosmic ray diffusion coefficients longitudinal 
and transverse relative to the magnetic field. 
For the simulations a Monte Carlo technique is 
used for the scattering of the particles in the up- 
stream, downstream and precursor media. The 
downstream spatial boundary required can be 
estimated from the solution of the convection- 
diffusion equation in a non-relativistic, large- 
angle scattering approximation in the downstream 
plasma which gives the chance of return to the 
shock, exp{—V2rb/d). The mean free path is 
calculated in the respective fluid rest frames (up- 



Figure 2: Here we show the acceleration time ra- 
tio (simulation time / theoretical time) versus the 
shock velocity for a representative shock incUna- 
tion of 75 degrees. We get similar results for other 
angles and as long as the shocks are oblique. The 
measurements are taken in the downstream side of 
the shock rest frame. Three runs have been per- 
formed. Starting from the bottom, the first Une is 
for a no precursor shock, the second with a pre- 
cursor width equal to 10 • A and the third with a 
precursor width equal to 100 • A. As one sees there 
is generally a reduction of the acceleration time as 
higher the speed of the shock, in any case, follow- 
ing the same trend of the relativistic shock acceler- 
ation simulations of past works. 

stream or downstream), assuming a momentum 
dependence to this mean free path for scattering 
along the field, related to the spatial diffusion co- 
efficient K, as we mentioned before. In the simu- 
lation the area of the precursor has a lower diffu- 
sion coefficient (where K|| >> kj^ is used (follow- 
ing Berezhko 1987) which physically reflects on an 
environment of a high density and large instabili- 
ties. The basic point in this work is that by chang- 
ing the diffusion conditions (by changing the val- 
ues of n) within the precursor of a width of many 
times the mean free path, one can study the scatter- 
ing of the particles which will physically reflect on 
the irregularity of the magnetic field. One impor- 
tant point to mention is that the adiabatic invari- 
ant is conserved within the precursor since in the 
de Hofmann-TeUer frame the adiabatic invariant is 
conserved anyway by crossing the shock from up- 
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stream to downstream and vice versa. Due to this 
condition, one can calculate a new pitch angle and 
velocity of the particle and the diffusion is continu- 
ing until the particle exits the precursor and conse- 
quently the shock region. Moreover, on the scatter- 
ing of the particle inside the precursor, we clarify 
two things: The mean free path A is dependent on 
the distance d inside the precursor in the following 
ways: 1) the diffusion coefficient k is a function of 
the position and 2) the mean free path is a function 
of energy. This means that as the particle moves 
in the precursor, its energy will change and hence 
its mean free path will change as well. The proba- 
bility that a particle will scatter after moving a dis- 
tance Ad, is Ad/ A. Comparing this ratio with a 
random number, defines as if the particle will con- 
tinue to scatter inside the precursor or exit this re- 
gion downstream or upstream. 
In figure 1 we show the shock velocity versus the 
obtained spectral index for three different simula- 
tion runs for a shock inclination equal to 75 de- 
grees. The first is for a shock with no precursor, 
the second is a shock with a precursor width equal 
to 10 • A and the third one with a a precursor width 
equal to 100 • A. One can see that the modified 
shocks create flatter spectra which at the 100 mean 
free paths width, is flatter then the sharp shock dis- 
continuity and the one with 10 mean free paths 
width. Also, we may see that all spectral indexes 
flatten as the shock velocity increases. In figure 
2 the acceleration rate versus the spectral index is 
shown for a shock incUnation of 75 degrees. Three 
runs have been performed. The first with no pre- 
cursor, the second with a precursor width equal to 
10 • A and the third with a precursor width equal 
to 100 • A. As one sees there is a reduction of the 
acceleration time. Furthermore, there is a trend for 
a reduction of the acceleration time as higher the 
speed of the shock, making it in any case an effi- 
cient cosmic ray accelerator. 

Conclusions 

We have preliminarily performed numerical sim- 
ulations, studying the behaviour of the accelerat- 
ing cosmic rays at highly obUque modified shocks, 
simulating the percursor with a length scale, L = 
k/V, corresponding to the diffusive scale of the 
energetic particles. It is clear that these modi- 



fied shocks result in flatter spectra compared to the 
pressumed sharp discontinuity shock fronts. This 
first finding is in accordance with theoretical pre- 
dictions, since the scattering inside the precursor 
confines high energy particles to further scattering, 
resulting in higher energies making the whole pro- 
cess more efficient. Moreover, one may see that 
the spectral index flattens as the shock velocity in- 
creases, leading to spectral index values of those 
compared to relativistic shocks. It also seems that 
even with modified shock fronts, there is an accel- 
eration speed up of the process as higher the speed 
of the modified shock (for speeds > 0.3c) and 
seem to be as efficient as the relativistic ones, pre- 
sented by [8], [9]. Additionally, this work has fur- 
ther implications connecting to the results of [10], 
where it is shown that even for the test-particle ap- 
proximation, highly oblique shocks can be very ef- 
ficient cosmic ray accelerators, achieving energies 
of lO^'^ eV. There is under way a detailed pre- 
sentation of this work. 
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